General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



ARL-TR-85-9 


cr- n i trf 

c.\ 

Copy No. 3 


COMPUTER SIMULATION OF THE EFFECTS OF A DISTRIBUTED 
ARRAY ANTENNA ON SYNTHETIC APERTURE RADAR IMAGES 
FINAL REPORT UNDER CONTRACT NAS9-17228 


James M. Estes 


APPLIED RESEARCH LABORATORIES 

THE UNIVERSITY OF TEXAS AT AUSTIN 
POST OFFICE BOX S029. AUSTIN, TEXAS 7S713-S02S 


9 April 1985 
Final Report 

1 September 1984 - 15 March 1965 

Approved for public release; 
distribution unlimited. 

Prepared for: 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
LYNDON B. JOHNSON SPACE CENTER 



(NASA-CB-17 Id86) CCflfUTfR SIHULATIOH Of THE N85-2C176 

EffECTS Of A DISIBlbl 1ED AfaHAl AN1EKNA Ch 
EiNTHE'lIC APEbTU bf biUf IflACiS final 

Report, 1 Eep. 1 5 fc 4 - 15 flat. 1585 (Texas Unclas 

Ucxv.) 28 p HC AC3/Hf A0 1 CSCL 171 G3/32 2 1740 


UNCLASSIFIED 


SECURITY C-ASSl FI CATION of THIS PACE (Whan Data Enfrma) 


REPORT DOCUMENTATION PAGE 

READ INSTRUCTIONS 
BEFORE COMPLETING FORM 

1. REPORT NUMBER 

2. GOVT ACCESSION NO. 

I. RECIPIENT'S CATALOG NUMBER 

4. TITLE (and Subtlt'a) 

COMPUTER SIMULATION OF THE EFFECTS OF A 
DISTRIBUTED ARRAY ANTENNA ON SYNTHETIC APERTURE 
RADAR IMAGES 

Final Report under Contract NAS9-17228 

S. TYPE OF REPORT • PERIOO COVERED 

final report 
1 Sep 84 - 15 Mar 85 

S. PERFORMING ORG. REPORT NUMBER 

ARL-TR-85-9 

7. Aj THORf a) 

James M. Estes 

S. CONTRACT OR GRANT HUMBERT*; 

NAS9-17228 

• PERFORMING ORGANIZATION NAME AND ADDRESS 

Applied Research Laboratories 
The University of Texas at Austin 
Austin, TX 78713-8029 

10. PROGRAM ELEMENT. PROJECT. TASK 
AREA • WORK UNIT NUMBERS 

11. CONTROLLING OFFICE NAME AND ADDRESS 

National Aeronautics and Space Administration 
Lyndon B. Jonnson Space Center 
Houston, Texas 77058 

12. REPORT DATE 

9 April 1985 

12. NUMBER OF PAGES 

27 

14 MONITORING AGENCY NAME 6 AOORESSf/f dlllerenl froai Cantrell Inf Olllce) 

IS. SECURITY CLASS, (ol I hie report) 

UNCLASSIFIED 

IS* OECL ASSI FICATION/ DOWNGRADING 
SCHEDULE 

<6. DISTRIBUTION STATEMENT (ol thl a Report) 

Approved for public release; distribution unlimited. 

>7. DISTRIBUTION STATEMENT (ol the abairacl aniarad in Block 30. II dltloranl from Report) 

IB. SUPPLEMENTARY NOTES 

synthetic aperture radar 
orbital SAR 
simulation 
SAR image analysis 

20 ABSTRACT (Continue on teveree aide II neceeeery end Identity by block number) 

The ARL:UT orbital SAR simulation has been upgraded to use three-dimensional 
antenna gain patterns. This report describes the modifications and presents 
quantitative image analyses of a simulation using antenna patterns generated 
from the modeling of a distributed array antenna. 


D!) iKw 1473 EDITION of I NOV 65 IS OBSOLETE UNCLASSIFIED 

SECURITY CLASSIF'CATION OF This PAOECWNm Data Knterod) 


















TABLE OF CONTENTS 


Page 


LIST OF FIGURES v 

LIST OF TABLES vil 

I. INTRODUCTION 1 

A. Background 1 

B. Distributed SAR Simulation 2 

II. IMPLEMENTATION 5 

A. OSS Modifications 5 

B. Image Analysis 5 

C. Curren': Simulation Model 6 

III. RESULTS 11 

A. 60° Nadir Angle Simulation 11 

B. 20° Nadir Angle Simulation 15 

C. 20° Nadir Angle Simulation, Short Array 15 

IV. CONCLUSIONS AND RECOMMENDATIONS 23 

REFERENCES 25 


PRECEDING P 


uj 


not ni vrn 


iii 



LIST OF FIGURES 


Figure Page 

1 Area of Terrain Model Imaged 7 

2 Complete Terrain Model 8 

3 SARCON Data 60° Nadir Angle Simulation 12 

4 SAR Images from 60° Nadir Angle Simulation 13 

5 SARCON Data 20^ Nadir Angle Simulation 16 

6 SAR Images from 20° Nadir Angle Simulation 17 

7 SARCON Data 20° Nadir Angle Simulation, 19 

Short Array 

8 SAR Images from 20° Nadir Angle Simulation, 20 

Short Array 


v 



LIST OF TABLES 


Table Page 

I PSL/NMSU Generated Antenna Patterns 3 

II Extracted Values from SAR Images for the 60° 

Nadir Angle Simulation 14 

III Extracted Values from SAR Images for the 20° 

Nadir Angle Simulation 18 

IV Extracted Values from SAR Images for the 20° 

Nadir Angle Simulation, Short Array 21 


vi i 



I. INTRODUCTION 


A. Background 

Under National Aeronautics and Space Administration Contracts 

NAS9- 15401 and NAS9-15217, Subtask 1, Applied Research Laboratories, The 

University of Texas at Austin (ARL:UT), developed a computer model of an 

orbiting synthetic aperture radar (SAR). This model enables the user to 

create a terrain model of point reflectors on an arbitrarily specifiable 

planet. The antenna is modeled by its gain pattern in azimuth and 

elevation and is fixed to the orbiter in accordance with the platform 

attitude and rates of change of attitude. The transmitter and receiver 

are specified by frequency, pulse repetition rate (PRF), impulse response 

characteristics, sampling rate, and range pulse compression codes. 

Orbital conditions are specified, as well as the area on the planet to be 

mapped and the geometry of the radar with respect to the area to be 

mapped. Synthetic in-phase and quadrature (IPQ) video data are then 

produced as the convolution of received scatterer echoes with the 

4 

transmitted waveform. The effects of antenna illumination, R 
attenuation, ambiguous ranging, and Doppler sampling are all taken into 
account. The IPQ video data are then processed through an SAR processor, 
and the resulting digital images are available for display and analysis. 
While the scope of the model is limited in some respects. It provides a 
powerful tool for certain parametric analyses (e.g., the effects of a 
parameter or part of an SAR may be investigated in an ideal system). 

The computer model developed and documented* under the 

aforementioned contracts is known as the orbital SAR simulation (OSS). 
Currently, the OSS consists of approximately fifty subroutines and 
programs and over 4500 lines of code written in FORTRAN IV and 
implemented on the ARL:UT CYBER 180/830. Since the completion of the OSS 
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version documented in Ref. 1, the OSS has been adapted to two subsequent 
2 3 

tasks/’- 3 

B. Distributed SAR Simulation 


This report under NAS9-17228 describes the effort to assess the 
effects of degradation in a distributed array antenna on SAR images. 
This study was performed in conjunction with the Physical Science 
Laboratory, New Mexico State University (PSL/NMSU). PSL/NMSU studied the 
effects of manufacturing tolerances of a dirtributed array antenna on its 
gain pattern. Amplitude and phase variations in the transmit/receive 
module were considered as were mechanical deformations. Then PSL/NMSU 
generated computer compatible tapes with three-dimensional gain patterns 
from arrays with various levels of degradations. 

The OSS was modified to use the antenna patterns generated by 
PSL/NMSU. A series of simulations representative of the shuttle imaging 
radar (SIR-C) were then executed. The antenna studies and the 
simulations were performed at C-band, 5.3 GHz. Two computer compatible 
tapes containing five antenna patterns each were received from PSL/NMSU; 
a summary of the data received is presented in Table I. The approach 
taken in this study was to set up a simulation and execute it using the 
ideal antenna pattern (tape DSAR2, file-1). The same simulation was then 
run with degraded antenna patterns and the resulting images were 
compared. 

In previous SAR studies ARL : UT adopted a set of "image quality" 
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parameters to characterize an SAR image. Most of these parameters 
(i.e., resolution, geometric fidelity, crispness, speckle, etc.) are a 
function of orbit geometry, radar hardware, or SAR processing parameters, 
and they would not be affected by variations in the antenna gain 
patterns. However, the extraction of backscatter levels (a 0 ) from the 
terrain is very likely to be affected by degradations in the antenna gain 
patterns. By placing fields with known values of a 0 in the terrain and 
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TABLE I 

PSL/NMSU GENERATED ANTENNA PATTERNS 
Data Files on Tape DSAR1 


Data File No. 

Description 

Gain (dB) 

1 

48 x 12 Array 22 dB Taylor elevation 
Taper - Ideal case 

35.29 

2 

48 x 12 Array 22 dB Taylor elevation 
Taper - parabolic bcw only 0.024 cm 

35.29 

3 

48 x 12 Array 22 dB Taylor elevation 
Taper - random T/R module variations 
at the element level 

Amplifier variations +1.76 dB , -3.0 d9 
Phase Variations ♦/- 10 

34.84 

4 

48 x 12 Array 22 dB Taylor elevation 
Taper - random T/R module variations 
(T/R module drives entire row) 

Amplifier variations +1.76 dB , -3.0 dB 
Phase Amplifier +/- 10 

35.11 

5 

48 x 12 Array 22 dB Taylor elevation 
Taper - parabolic bow (0.024 cm) and 
random T/R module variations at the 
element level 

Amplifier variations +1.76 dB, -3.0 dB 
Phase variations +/- 10 

34.81 


Data Files on Tupe DSAR2 



288 x 12 Element DSAR Array 
22 dB Taylor Elevation Taper 
Array dimensions 12.4 x 0.54 m 


Data File No. 

Description 

Gain (dB) 

1 

Ideal case - no mechanical deformations 
No T/R module degradations 

43.98 

2 

Bow only - 0.175 cm oarabollc bow 
No T/R module degradations 

43.96 

3 

Weight only - no mechanical deformations 
T/R module variations: amplitude +1.76, 
-3 dB, phase +/- 10 
288 x 12 T/R modules (at element level) 

43.63 

4 

Both - 0.175 cm parabolic bow 

T/R module variations: amplitude +1.76, 

■•3 dB, phase +/- 10 

288 x 12 T/R modules (at element level) 

43.59 

5 

Both 2 - 0.175 cm parabolic bow 

T/R mouule variations: amplitude +1.76, 

-3 dB, phase +/- 10 

16 x 12 T/R modules (each module feeds a 

row, In a sub-panel, 18 elements long) 

43.53 




then extracting the val'es for o° from the images, the effects of the 
antenna degradations on the measurement of o° can be assessed. It is 
important at this point to note that the shape of the antenna gain 
pattern affects the measuremjnt of o° within an image as a function of 
both range and azimuth. If the exact pattern is known, however, then the 
values for o° may be properly adjusted (assuming there is no asymmetrical 
structure). We are concerned here with the degradations which cause the 
antenna pattern to deviate from what was assumed. In other words, we are 
interested in the effects on an SAR image due to the antenna degradations 
after the pattern has been measu.ed. 
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II. IMPLEMENTATION 


A. OSS Modifications 


The basic change required to the OSC was a modification of the 
program which uses the antenna pattern (program ECHO). The original 
simulation used elevation and azimuth principal plane cuts to model the 
antenna. Program ECHO and subroutine ANTWT were modified to use the 
three-dimensional antenna patterns generated by PSL/NMSU (redesignated 
PECHO and PANTWT). The antenna patterns were comprised of 181 conical 
cuts (361 samples per cut, 1° in 4> per sample) of 0.5° or 0.25° in 0 per 
cut. Since the 3 dB beamwidth in azimuth of the full array was on the 
order of 0.25°, an interpolation scheme was also Included in program 
PECHO. For efficiency, the original version of the OSS did not use an 
interpolation scheme; instead the desired antenna pattern was 

sufficiently sampled so that interpolation was not necessary. However, 
In the present case, sufficient sampling of the antenna pattern would 
have generated more data than could fit into memory. A quadratic fit 
through three points was used for interpolation. 

A program (ANTPSL) was also written to generate the file for use by 
program PECHO. ANTPSL takes the antenna patterns generated by PSL/NMSU 
on 9-track computer compatible tapes and decodes them into a binary file. 
More operator interaction is required in the modified antenna simulation 
than in the original model, however, because the modified antenna 
simulation was not fully integrated into the OSS. 

B. Image Analysis 

The current version of the OSS provides no software for the 
extraction of Image statistics needed for Image quality analyses. The 
only post-processing routines available (POST) generated a tape for off- 
line image display on an ARL:UT graphics system, which is no longer 
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available. A version of program POST (PPOST) was generated which creates 
a computer tape of the output SAR Image for use on an HP 21MX 
minicomputer with a urlnnel color graphics display (we used only gray 
scale and an attached Tektronix hard copy unit). 

Another version of program POST (PPOSTST) was created which 
generates a binary file of the SAR Image. The structure of the file is 

compatible with that generated by an SAR processor which uses real radar 
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data. Programs were available from previous work which allowed 

extraction of statistics (max, min, mean, and standard deviation) from 

the SAR Images. For this study, these programs were used to print the 

SAR image pixel values from which the locations of the fields in the 

Images could be deduced. Next, the portions of the Images containing the 

fields were statistically analyzed, and normalized cross -section 

measurements were obtained. 

C. Current Simulation Model 


For the simulations performed, the planet modeled was earth. The 
orbit parameters were chosen to correspond as closely as possible with 
the shuttle's, 250 km altitude, near-polar orbit. The frequency was 
C-band (5.3 GHz), the ground range resolution was 25 m with a sample 
ratio of 2 (i.e., sample interval 1/2 of resolution), and the overlay 
ratio was 2 (two "looks"). Two simulations were performed at a 20° nadir 
angle and one at a 60° nadir angle. 

The area of the terrain model to be imaged is shown in Fig. 1. For 
the 20° nadir angle simulation (using the full array) the ambiguity 
locations were computed (see p. 4/60 of Ref. 1 for a discussion of 
ambiguity calculations); fields were then placed at these locations so 
that the ambiguities would be introduced into the images. The antenna 
pattern and SAR parameters act to suppress these ambiguities. The 
complete terrain model, used for all three simulations, is shown in 
Fig. 2. The ambiguity locations for the 60° nadir angle simulation are 


6 


ORIG'VA’l. P' ■ 

OF POCiu . *1.1 < i 



FIGURE 1 

AREA OF TERRAIN MODEL IMAGED 
ALL FIELDS CONTAIN 4 SCATTERERS PER RESOLUTION CELL 
(25 m x 2? m), WHICH ARE RANDOMLY DISTRIBUTED 


ARL:UT 
AS-85-407 
JME GA 
4-10 85 


7 



AMBIGUITY 

FIELDS 


PATCH 

CENTER 



ARL UT 
AS-86-408 
JMC • GA 
4-10 65, 


8 


different than for the 20° r.adir angle simulation. However, the 
unambiguous sector imaged was small, and some of the fields in ambiguity 
locations fall in the azimuth sidelobes of the antenna pattern. 
Therefore, the same terrain model could be used for all the simulations. 
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III. RESULTS 


A. 60° Nadir Angle Simulatio n 

The parameters used for this simulation are shown in Fig. 3. The 
simulation was exercised for the antenna patterns from tape DSAR2, 
files 1 and 5 (see fable I). The SAR images produced by the OSS are 
shown In Fig. 4. The normalized cross-section values from the imaged 
fields are shown in Table II for the two antenna patterns. Two areas, A1 
and A2, are outlined in the SAR images of Fig. 4. These are areas where 
no targets were placed; the values obtained for A1 and A2 represent 
ambiguously imaged fields. 

The variation in the computed scattering levels for imaged fields 
within one image is a function of the antenna gain pattern shape and the 
portion of the radat at the start of sampling. The two-way, 3 d5 width 
of the distributed array is approximately 0.2°; consequently, at a slant 
range of 530 km, only 1.8 km of the terrain is illuminated. The 
synthetic array length was 0.7 km. In a real system these variations 
could be accounted for and need not concern us here. The differences in 
the values obtained for the ideal and degraded patterns do not take into 
account the difference in the gain of the two patterns indicated in 
Table I. For full array patterns (DSAR2) the difference in gain between 
the two patterns used was 0.45 dB; this translates to a 0.9 dB overall 
reduction in the degraded pattern results. If a calibration were not 
performed after deployment of the antenna, this could be significant. 


The values obtained for sections A1 and A2 give some idea of the 
ambiguity level to be expected. It is important to note that the terrain 
illuminated by the antenna sidelobes was not filled with targets, and 
therefore the values for A1 and A2 do not represent the integrated 
sidelobe ratio. However, the fields in the terrain model do fall in the 
first few major sidelobes, and it is reasonable to expect that the values 
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ORIGINAL PAG- IS 
OF POOR QUALITY 


PLANET NAME ** 
ECCENTRICITY ** 
ROTATIONAL RATE (0E0/8) 


PLANET SPECIFICATION • *•#** * * ##* 

EARTH EQUATORIAL RADIUS (KM) *• . 63781 670000E 
. 018201 79996E-01 GRAVITATIONAL CONSTANT (KM3/SEC2) ** 39860 lOOOOOE 
. 41 78074 9999E-02 TIME OF PRIME MERIDIAN PASSAGE (S) **0. 


ORBIT SPECIF.' CATION 


ORBIT I. D. ** 

ECCENTRICITY *# 

LONO OF ASCEND I NO NODE (DEO) ** 
TIME OF PERIOEE PASSAGE (S) *• 

ORBITER INITIALIZATION TIME(S) ** 


SHUTTLE SEMI -MAJOR AXIS (KM) ** 

. 20000000000E-02 INCLINATION (DEG) ** 

0. ARGUMENT OF PERIGEE (DEG) •* 

0. ROTATIONAL RATE (DEO/S) ** 

1 OOOOOOOOOOE-03 


. 66280000000E 
. 1 0800000000E 


0 . 


. 67037940870E 


RADAR ID. •• 

RECEIVER/TRANSMITTER BW (MhZ) ** 
SIONAL-TO-NOISE RATIO (DB) ** 
SAMPLE LENOTH OF RANOE CORRELATION 
BINARY PHASE CODE •• 

GROUND RANOE RESOLUTION (M) •• 

RANOE SAMPLINO RATIO (M) ** 

RANOE IMPULSE RESPONSE FUNCTION •« 
PATCH-TO-PATCH OFFSET. RNO (M) ** 

RANOE SWATH WIDTH (KM) ** 

MAP 8TART LATITUDE (DEO) •• 


RADAR SPEC 
SIR-C 

. 66609 1 998 13E 4-01 
. lOOOOOOOOOOE^OS 
1 

B 

. 29000000000E+02 
. 20000000000E+0 1 
C0SINE**2 


0. 


. 40000000000E+01 
. 49000000000E+02 


IF ICATION ********4 
OPERATING WAVELENGTH (M) ** 

RANCE TIME-BANDWIDTH PRODUCT *# 

A/D 8AMPLE RATE (MHZ) •« 

SAMPLE LENOTH ACROSS PHASE INTERVAL 
BINARY PHASE CODE SEQUENCE •• 
AZIMUTH RESOLUTION (M) ** 

AZIMUTH SAMPLINO RATIO (M) ** 

APERTURE WEIGHT FUNCTION ** 
PATCH-TO-PATCH OFFSET. AZ (M) ## 

NO OF PATCHES ** 


. 96600000000E 
. 10000000000E 
. 13321031 163E 

RANDOM 

. 2 9000000000E 
. 20 OOOOOOOOOE 
TAYLOR 


0. 


ANTENNA SPECIFICATION 


ANTENNA ID. •» 

BORESIOHT SQUINT AT TO (DEO) ** 
ELEVATION ANGULAR COVERAGE (DEO) 
PHASE CENTER. BODY AXIS Y (M) •* 

COORD SYS. BODY AXIS ROLL. (DEO) 
COORD SYS. BODY AXIS YAW. (DEO) 
PLAT PITCH RATE (DEO/S) ** 


6-PANEL-C BORESIOHT NADIR AT TO (DEO) *# 

. 90000000 OOOE +02 AZIMUTH ANOULAR COVERAOE (DEO) • 
** . 62000000000E+01 PHASE CENTER. BODY AXIS X (M> ** 
0. PHASE CENTER. BODY AXIS Z (M) ** 

** 0. COORD SYS. BODY AXIS PITCH 

»* . 90000000 ")00E402 PLAT ROLL RATE (DEO/S) ** 


. 60000000000E- 
. 30000000000E ; 

° 

A 


(DEO) ** -. 30000000000E- 

0. 


0. 


PLAT YAW RATE (PEG/S) ** 


TERRAIN I. D. •• 

NO OF FIELDS ** 

X-AXIS COVERAOE (KM) ** 
TERRAIN CENTER. R (KM) ** 
TERRAIN CENTER. LONO (DEO) ** 


TERRAIN SPECIFICATION * * * * 
ANT-EVAL NO OF DISCRETES *# 

13 TOTAL NO OF SCATTERERE ** 

. 3 1 OOOOOOOOOE 4-02 Y-AXIS COVERAOE (KM) ** 

. 63679 192370E 4-04 TERRAIN CENTER. LAT (DEC) ** 
. 34463288002E4-03 


• »•••• SYNT 
SYNTHETIC ARRAY NO ** 

ARRAY LENOTH (M) #* 

ARRAY FORMATION TIME (MS) ** 

NO OF PULSES •* 

NO OF RANGE SAMPLES ** 

PATCH CENTER RANOE SAMPLE NO. ** 
SLANT RANGE RESOLUTION (M) »• 

START RANOE (KM) #• 

SQUINT ANOLE PATCH CENTER (DEO) *« 
LOS AZIMUTH AT PATCH CENTER (DEC) ** 
PATCH CENTER. R (KM) ** 

PATCH CENTER. LAT (DEO) ** 

PATCH CENTER. LONO (DEO) ** 


H E T I C A R R 
1 

. 68480971 927E+03 
. 881B3493374E4-02 
138 
319 
160 

. 229092022 13E4-02 
. 93479779887E4-03 
. 9004U60378E+02 
-. 1 1 1367224 16E4-03 
. 63679 192370E4-04 
49000000000E4-02 
34463288002E4-03 


AY PARAMETERS ***** 
TRANSMISSION START TIME (S) *• 

ARRAY INCLINATION (DEG) •• 

PLATFORM VELOCITY (KM/S) »* 

PRF (HZ) ** 

NO. OF AZIMUTH FILTERS ** 

SLANT RANGE SWATH WIDTH (KM) •• 

SLANT RANGE SAMPLE INTERVAL (M) ** 

RANCE PATCH CENTER (KM) ** 

NADIR ANGLE PATCH CENTER (DEG) •• 

LOS INCIDENCE AT PATCH CENTER (DEG)** 
ORBITER MASS CENTER. R (KM) •• 
ORBITER MASS CENTER. LAT (DEO) ** 
ORBITER MASS CENTER, LONO (DEG) •« 


7< 

. 1 1 900000000E 
. 49000000000E 1 


. 68663934395E". 
. 82991 9911 60E- 
. 77696899903E-; 
. 1 9637937780E- 

. 3600826633 IE 
. 1 1292601 106E- 
93654696249E- 
. 60000986690E- 
64 1 89433038E- 
66 1881 09998E' 

. 43347470309E 
. 33927066739E- 


FIGURE 3 

SARCON DATA 60° NADIR ANGLE SIMULATION 
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TABLE II 


EXTRACTED VALUES FROM SAR IMAGES 
FOR THE 60° NADIR ANGLE SIMULATION 


FIELD 


IDEAL PATTERN DEGRADED PATTERN 

DSAR2/Fi ? e 1 DSAR2/ Fi 1 e 5 

(dD) CdB) 


- 3 dB 

0.2 

0.16 

- 6 dB 

- 4.31 

- 4.35 

- 9 dB 

- 5.33 

- 5.36 

-12 dB 

- 9.61 

- 9.61 

-15 dB 

-11.41 

-11.42 

+ 3 dB 

+ 5.96 

+ 5.92 

A1 

-12.3 

-12.0 

A2 

-16.0 

-15.6 

Patch Center 

- 1.79 

- 1.79 
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for A1 and A2 are within 2-3 dB of the total ambiguity level. From these 
data we can infer a value for dark target contrast (ratio of the mean 
level of an area of very low reflectance to the mean level of the entire 
image) of approximately -14 dp. It appears that most of the energy is 
from a field ambiguously imaged in the main lobe. 

B . 20° Nadir Angle S imulation 

Tne parameters used for this simulation are shown in Fig. 5. The 
simulation was performed using files 1 and 5 from 0SAR2. The SAR images 
produced are shown in Fig. 6. Table III is a list of the extracted 
values from the images. As pointed out in Section III. A, the variation 
in reflectivity level, in azimuth, across the image is due to the shape 
of the antenna gain pattern. 

For this simulation, the values for ambiguous areas A1 and A2 are 
much lower than in the previous simu’ation. The pulse repetition 
frequency (PRF) for this simulation unambiguously images 0.5° at the main 
beam, instead of the 0.3° imaged in the previous simulation. 

C. 20° Nadir Angle Simulation, Short Array 

The geometry of this simulation was identical to the previous 20° 
nadir angle case. The antenna patterns used, however, were from tape 
DSAR1 (see Table I). These patterns were created from one panel (out of 
six) of the distributed array antenna and have a corresponding beamwidth 
six times that of the full array in azimuth (elevation dimensions are 
identical). To avoid ambiguity problems, a higher PRF was used. 

The parameters used for this simulation are shown in Fig. 7. The 
SAR images produced using the antenna patterns from DSAR1 , files 1, 4, 
and 5, are shown in Fig. 8. The corresponding values for normalized 
cross-section are presented in Table IV. 
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«**«*«,*« PLANET SPECIFICATION **«##*##• 

PLANET NAME •• EARTH EQUATORIAL RADIUS (KM) ** . 63781670000E 

ECCENTRICITY ** . 81820179996E-01 GRAVITATIONAL CONSTANT (KM3/SEC2) ** . 39860 lOOOOOE 

ROTATIONAL RATE (DEO/8) ** . 41780749999E-02 TIME OF PRIME MERIDIAN PA8SA0E (S) **0. 

• •»•••••• ORBIT SPECIFICATION ********** 

ORBIT I.D. •• SHUTTLE SEMI-MAJOR AXIS (KM) ** . 66280000000E 

ECCFNTRICITY ** . 20000000000E-02 INCLINATION (DEO) ** . 10800000000E 

LONG OF ASCENDINO NODE (DEO) •• 0. ARGUMENT OF PERIGEE (DEO) ** 0. 

TIME OF PERIOEE PASSAGE (S) ** 0. ROTATIONAL RATE (DEO/S) •• . 67037 940870E 

ORB I TER INITIALIZATION TIME(8> ** -. 10000000000E-03 


• ••••••*« RADAR SPECIFICATION ********** 

RADAR I.D. *# SIR-C OPERATING WAVELENOTH (M) ** . 96600000000E 

RECEIVER/TRANSMITTER B W (MHZ) ** . 16863430288E+02 RANGE TIME-BANDWIDTH PRODUCT ** . 10000000000E 

S I GNAL-TO-NO I SE RATIO (DB) ** . 10000000000E+03 A/D SAMPLE RATE (MHZ) ** . 33726860976E 

SAMPLE LENOTH OF RANGE CORRELATION ** 1 SAMPLE LENGTH ACROSS PHASE INTERVAL * 

BINARY PHASE CODE ♦* B BINARY PHASE CODE SEQUENCE •* RANDOM 

GROUND RANGE RESOLUTION (M) •• 29000000000E+02 AZIMUTH RESOLUTION (M) *« . 29000000000E 

RANGE SAMPLING RATIO (M> •• . 20000000000E+01 AZIMUTH SAMPLING RATIO (M) •• . 20000000000E 

RANGE IMPULSE RESPONSE FUNCTION •• C0SINE#*2 APERTURE WEIOHT FUNCTION ** TAYLOR 

PATCH-TO-PATCH OFFSET. RNO (M) ** 0. PATCH-TO-PATCH OFFSET. AZ (M) ** O. 

RANOE SWATH WIDTH (KM) •• ' . 40000000000E+01 NO OF PATCHES •• 


MAP START LATITUDE (DEO) •• . 49000000000E+02 


• ••#«»•• ANfENNA SPECIFICATION *•#•••••• 

ANTENNA I.D. ** 6-PANEL-C BORESIGHT NADIR AT TO (DEG) •• . 20000000000E 

BORESIOHT SQUINT AT TO (DEO) ** . 90000000000E+02 AZIMUTH ANOULAR COVERAGE (DEO) •• . 90000000000F 

ELEVATION ANOULAR COVERAGE (DEG) ** . 62000000000E+01 PHASE CENTER. BODY AXIS X (M) ** 34000000000E 

PHASE CENTER. BODY AXIS Y (M) #* 0. PHASE CENTER. BODY AXIS Z (M) *• 0. 

COORD SYS. BODY AXIS ROLL, (DEG) ** 0. COORD SYS. BODY AXIS PITCH, (DEO) •• -. 70000000000E 

COORD SYS. BODY AXIS YAW. (DEO) •• 90000000000E+02 PLAT ROLL RATE (DEO/S) »• 0. 

PLAT PITCH RATE (DEO/S) ** 0. PLAT YAW RATE (DEG/S) ** 0. 

• ••••*«• TERRAIN SPECIFICATION • •»•••••* 

TERRAIN I. D. ** ANT-EVAL NO OF DISCRETES ** 

NO OF FIELDS ** 13 TOTAL NO OF SCATTERERS ** 7 

X-AXIS COVERAGE (KM) ** . 31000000000E+02 Y-AXIS COVERACE (KM) *# . 1 1 900000000E 

TERRAIN CENTER, R (KM) ** . 636791 92370E+04 TERRAIN CENTER. LAT (DEO; #* . 49000000000E 

TERRAIN CENTER, LONG (DEO) ** . 33937463699E+03 


SYNTHETIC ARRAY PARAMETERS *«•••• 

SYNTHETIC ARRAY NO ** 1 TRANSMISSION START TIME (S) ** . 7080761 7297E 
ARRAY LENOTH (M) •• . 34254838328E+03 ARRAY INCLINATION (DEG) •• . 849984 18492E 
ARRAY FORMATION TIME (MS) ## . 44 1 1210206 7E+02 PLATFORM VELOCITY (KM/S) *• . 7769409731 9E 
NO OF PULSES *# 119 PRF (HZ) *# . 26090123991E/ 
NO OF RANOE SAMPLES *# 319 NO. OF AZIMUTH FILTERS *# 

PATCH CENTER RANOE SAMPLE NO. •• 160 SLANT RANGE SWATH WIDTH (KM) •• . 14221779309E 
SLANT RANOE RESOLUTION (M> ## . 88888390377E+01 SLANT RANOE SAMPLE INTERVAL (M) . 44444179188E- 
START RANOE (KM) •• . 26767702799E+03 RANGE PATCH CENTER (KM) *• . 2A838368998E 
SQUINT ANGLE PATCH CENTER (DEO) •« . 90104737938E+02 NADIR ANGLE PATCH CENTER (DEO) •* . 2000 1321 723E 
LOS AZIMUTH AT PATCH CEfiTER (DEO) •• -. 1 1 489683437E+03 LOS INCIDENCE AT PATCH CENTER (DEO)** . 20827364886E 
PATCH CENTER. R (KM) •• . 63679 192370E+04 ORDITER MASS CENTER. R (KM) »• . 66190939439E 
PATCH CENTER. LAT (DEG) •• 45000000000E + 02 ORB ITER MASS CENTER. LAT (DEO) •• . 44647387989E 
PATCH CENTER. LONG (DEG) •• . 33937463699E+03 ORBITER MASS CENTER. LONG (DEG) *• 338321 43834E- 


FIGURE 5 

SARCON DATA 20° NADIR ANGLE SIMULATION 
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TABLE III 


EXTRACTED VALUES FROM SAR IMAGES 
FOR THE 20° NADIR ANGLE SIMULATION 


FIELD 

IDEAL PATTERN 
DSAR2/File 1 
(dB) 

DEGRADED PATTERN 
DSAR2/ File 5 
(dB) 

- 3 

dB 

- 0.02 

0.41 

- 6 

dB 

- 7.13 

- 6.94 

- 9 

dB 

- 5.63 

- 5.96 

-12 

dB 

-12.63 

-12.35 

-15 

dB 

-11.54 

-11.51 

+ 3 

dB 



A1 


-27.6 

-27.8 

A2 


-51.8 

-50.7 

Patch Center 

- 2.02 

- 2.01 
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fci.NET NAME *• 
ENTRICITY •• 

AT I O' ML RATE (DEO/S) 


ORIGINAL PAGE IS 
OF POOR QUALITY 

PLANET SPECIFICATION *••****•* 

EARTH EQUATORIAL RADIUS (KM) ** * . 637B1670000E+04 

. 010201 79996E-01 GRAVITATIONAL CONSTANT (KM3/SEC2) *• . 39860 100000E+06 
. 41780743999E-02 TIME OF PRIME MERIDIAN PA88A0E (S) **0. 


ORBIT SPECIFICATION 


IT I. D. ** 

ENTRICITY ** 

(0 OF A8CENDIN0 NODE (DEO) •• 

E OF PERIOEE PASSAGE (8) ** 

ITER INITIALIZATION TIME(3) •» 


8HUTTLE SEMI -MAJOR AXIS (KM) *• 

. 20000000000E-02 INCLINATION (DEO) •• 

0. ARGUMENT OF PERIOEE (DEG) ** 

0. ROTATIONAL RATE (DEO/8) #* 

1 0000000000E-03 


66280000000E+04 
. 1 0800000000E+03 
). 

. 67037940870E-0 1 


•AR I D. •* 

EIVER/TRAN8MITTER BW (MHZ) *• 
'NAL-TO-NO I BE RATIO (DB> »« 

IPLE LENOTH OF RANOE CORRELATION 
IIARY PHASE CODE *# 

JUNO RANOE RESOLUTION (M) ** 

JOE 8AMPLIN0 RATIO (M) •• 

iGE IMPULSE RESPONSE FUNCTION •• 
CH-TO-PATCH OFFSET, RNO (M) *# 

tOE SWATH WIDTH (KM) ** 

START LATITUDE (DEO) •• 


ADAR SPEC 
SIR— C 

. 1 68634 30288E +02 
. i 0000000000£+03 

•• 1 
B 

. 29000000000E+02 
. 20000000000E+0 1 
COS I NE#*2 
0. 

. 40000000000E+01 
. '♦3000000000C+02 


IF ICATION ******♦•« 
OPERATING WAVELENGTH (M> •• 

RANGE TIME-BANDWIDTH PRODUCT ** 

A/D 8AMPLE RATE (MHZ) ** 

SAMPLE LENOTH ACROSS PHASE INTERVAL 
BINARY PHASE CODE SEQUENCE •• 
AZIMUTH RESOLUTION (M) ** 

AZIMUTH 8AMPLIN0 RATIO (M> ** 

APERTURE WEIOHT FUNCTION *• 
PATCH-TO-PATCH OFFSET, AZ (M) #* 

NO OF PATCHES ** 


. 36600000000E-0 I 
. 10000000000E+01 
. 33726860976E+02 
1 

RANDOM 

. 2 90000 0000 OE >02 
. 20000000000E+01 
TAYLOR 


0. 


ANTENNA 


SPECIFICATION 


ENNA I. D. •• 
pESIOHl SQUINT AT TO (DEO) ** 
frVATION ANOULAR COVERAOE (DEO) 
fiSE CENTER, BODY AXIS Y (M) •• 

IRD SYS, BODY AXIS ROLL. (DEO) 
)RD SYS, BODY AXIS YAW, (DEO) « 
*T PITCH RATE (DEO/8) ** 


6-PANEL-C BORESIGHT NADIR AT TO (DEO) •• 

90000000000E+02 AZIMUTH ANOULAR COVERAOE (DEG) ♦♦ 
. 90000000000E+0 1 PHASE CENTER. BODY AXIS X (M) ** 
0. PHASE CENTER. BODY AXIS Z (M) •* 

0. COORD SYS. BODY AXI8 PITCH. 

. 90000000000E+02 PLAT ROLL RATE (DEO/S) •* 

0. PLAT YAW RATE (DEO/S) ** 


. 20000000000E+02 
. 1 9000000000E+01 
0. 

0. 

(DEO) ** 70000000000E+02 

0. 

0 . 


TERRAIN SPECIFICATION 


DRAIN I. D. •• 

OF FIELDS ** 
tXIS COVERAOE (KM) •* 

I (RAIN CENTER, R (KM) •* 
IRAIN CENTER, LONO (DEO) •• 


ANT-EVAL 

13 

. 3 1 000000000E+02 
. 63679 192370E4-04 
. 33937463695E+03 


NO OF DISCRETES ** 

TOTAL NO OF SCATTERERS ## 
Y-AXIS COVERAOE (KM) •• 
TERRAIN CENTER, LAT (DEO) ## 


71 

7099 

. 1 1 900000000E+03 
. 4 3000 00000 0E+02 


• ••••• SYNT 
JTHETIC ARRAY NO •* 

JAY LENOTH (M) ## 

1AY FORMATION TIME (MS) #* 

OF PULSES ** 

OF RANOE SAMPLES ** 

'CH CENTER RANOE SAMPLE NO. ** 

\NT RANOE RESOLUTION (M) ** 

1RT RANOE (KM) *• 

1INT ANGLE PATCH CENTER (DEO) *• 

5 AZIMUTH AT PATCH CENTER (DEO) ** 
rCH CENTER, R (KM) *# 
rCH CENTER, LAT (DEO) ** 

'CH CENTER, LONO (DEO) •• 


H E T I C A R R 
1 

. 34294838328E+03 
. 441121 02067E+02 
333 
319 
160 

. 88888390377E+01 
. 26767702799E+03 
. 90 1 04737338E+02 
-. 1 1489683437E+03 
. 63679 192370E +04 
. 43000000000E+02 
. 33937463693E+03 


AY PARAMETERS ***** 
TRANSMISSION START TIME (S) ** 

ARRAY INCLINATION (DEG) •* 

PLATFORM VELOCITY (KM/S) *• 

PRF (HZ) ** 

NO. OF AZIMUTH FILTERS ** 

SLANT RANGE SWATH WIDTH (KM) ** 

SLANT RANGE SAMPLE INTERVAL (M) ** 

RANOE PATCH CENTER (KM) ** 

NADIR ANGLE PATCH CENTER (DEO) ** 

LOS INCIDENCE AT PATCH CENTER (DEO)** 
ORB I TER MASS CENTER. R (KM) •• 

ORB I TER MASS CENTER. LAT (DEO) •• 

ORB I TER MASS CENTER. LONO (DEO) •• 


. 7080761 7297E+03 
. 849384 18492E-01 
. 7769409731 5E+01 
. 799493384 13E+04 
200 

14221773309E+01 
. 444441 79 188E+01 
26830368998E+O3 
. 20001 32 1723E+02 
. 20B27364886E+02 
. 66190939439E-* , 04 
. 44647387989E+02 
33832 143834E+03 


FIGURE 7 

SARCON DATA 20° NADIR ANGLE SIMULATION, SHORT ARRAY 
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TABLE IV 


F I ELD 


- 3 dB 

- 6 dB 

- 9 c'B 
-12 dB 
-15 dB 
+ 3 dB 

A1 

A? 

Patch Center 


EXTRACTED VALUES FROM SAR IMAGES 
FOR THE 20° NADIR ANGLE SIMULATION, SHORT ARRAY 


IDEAL PATTERN DEGRADED PATTERN DEGRADED PATTERN 
DSAR2/F1 1 e 1 DSARl/FIle 4 DSARl/Fi le 5 

(dB) [dB] (dB) 


+ 1.64 

+ 1.70 

+ 1.64 

- 1.68 

- 1.68 

- 1.70 

- 4.88 

- 4.86 

- 4.88 

- 7.53 

- 7.52 

- 7.53 

- 9.91 

- 9.90 

- 9.91 


+ 7.82 

+ 7.87 

-20.4 

-20.4 

-20.4 

-16.8 

-16.8 

-16.8 

+ 0.05 

+ 0.05 

+ 0.05 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


The OSS has been modified to use three-dimensional anLenna patterns 
(specified as conical cuts). The simulation analyzed in this report was 
exercised on three configurations with seven complete simulations. A 
copy of the OSS modifications and programs used to analyze the images has 
been sent to NASA/JSC on computer compatible tape, along with another 
copy of the original OSS. 

The results summarized in Tables 1 1 - IV show very little effect due 
to distortions in the shape of the antenna patterns. The effect on the 
images due to pattern distortions is on the order of a few tenths of a 
decibel. Additionally there is a reduction in overall gain of up to 
0.9 dB for the degraded patterns (see Table I for the one-way gain 
reductions). A comparison of the antenna patterns received from PSL/NMSU 
shows virtually no perceptible main beam distortion or pointing error. 

The results of this study show that, for the electrical variations 
in T/R modules on the order of +1.8/-3 dB and ±10° phase, the effects on 
the SAR images are negligible except for a gain factor which would 
presumably be measured after antenna construction and before deployment. 
Implicit in the previous statement is the assumption that the electrical 
variations a^e the result of manufacturing tolerances and are not time 
dependent. It is appropriate to note here that the antenna patterns were 
modeled as oain versus angle, and that if there are any phase distortions 
with angle across the main beam, then several parameters in an image 
including resolution and SAR beamforming could be affected. The 
mechanical deformations which were modeled, however, represent 
symmetrical variations in the manufacture of the antenna and may not be 
representative of errors in the deployment of the antenna in space. 
Asymmetrical mechanical deformations may be more serious than what have 
been model ad to date. 


23 


•,u.u:OiN<3 TA u: 


'! \NK N ■ • : 





Effects on tK SAR image not fuliy modeled in this study include the 
effects of amoiguities. The PRFs proposed for the SIR-C (1200- 
1900 pulses/s) will net allow unambiguous imaging o' the main lobe of the 
antenna pattern. In the 60° nadir simulation performed in this study, 
the PRF was 1570 Hz, which corresoonds to an unambiguous imaging coverage 
of 0.3°. The one-way 5 d3 width of the main lobe of the full array is 
approximately 0.3°, so that for this set-up ambiguo-s returns on the 
order of -12 dB would be observed. T he other mechanism for ambiguous 
imaging is the reception of scattered energy in the sidelobes of the 
antenna. Without shading, the azimuth sidelobes are on the order of 
-13 dB from the peak. Since ambiguous imaging through azimuth sidelobes 
occurs from both transmission and reception through the sidelobes, the 
received energy will be 26 dB below the peak. In the case of range 
ambiguities which occur by eclipsing of the received pulses, however, 
energy transmitted on several pulses is received through the rar.ge 
sidelobes of each pulse, thus effectively changing the range sidelobe 
pattern by up to 10 dB. The range sidelobes of the patterns generated by 
PSL/NMSU were 22-24 dB below the peak (22 dB Taylor weighting was used). 
For the single panel case generated by PSL/NMSU, with T/R modules driving 
rows of 48 elements, the range sidelobes were only 19 dB down from the 
peak. Depending on the configuration, these may be unacceptably high. 

Future studies could take two paths. One would be to more closely 
mode 1 the proposed SIR-C and investigate the effects of ambiguities and 
mechanical deformations of the antenna. The second line of research 
would be to investigate the possible effects of mechanical deformations 
due to deployment, thermal gradients, and other time varying factors on 
the SAR images. Since measurement of the antenna gain pattern afeer 
oepl^yment does not seem possible, these time varying tactors can create 
irrecoverable errors in the SAR images. 
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